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In recent articles [Mishima et ai, Phys. Rev. A, 66, 033401 (2002); Chao, Phys.Rev. A, 72, 
■ 053414 (2005)] it was proposed to use the residue theorem for the exact calculation of the tran- 

sition amplitude describing strong-field ionization of atomic systems within Keldysh theory. This 
should avoid the necessity to apply the method of steepest descent (saddle-point approximation). 
(N . Comparing the results of both approaches for atomic hydrogen a difference by a factor of 2 was 

found for the Is, and an even more drastic deviation for the 2s state. Thus it was concluded that 
the use of the saddle-point approximation is problematic. In this work the deviations are explained 
and it is shown that the previous conclusion is based on an unjustified neglect of an important 
' contribution occurring in the application of the residue theorem. Furthermore, the applicability of 

_ the method of steepest descent for the ionization of Rydberg states is discussed and an improvement 

of the standard result is suggested that successfully removes the otherwise drastic failure for large 
principal quantum numbers. 
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I. INTRODUCTION 

The ionization process in atomic systems exposed to 
a strong laser field has attracted great interest during 
the past decades and its proper modeling remains a chal- 
lenge to theory [l[ . Among the numerous approximations 
developed to treat this problem Keldysh theory Q pos- 
sesses a prominent role. This theory was proposed by 
Keldysh 40 years ago and describes the ionization pro- 
cess as a transition between an initial electronic bound 
and a Volkov continuum state (adopting the length gauge 
formulation). Besides the main approximation of the 
Keldysh theory, the neglect of the interaction of the es- 
caping electron with the long-range Coulomb potential 
(in the case of a neutral atom), two additional simpli- 
fications were made in 0]: (i) the method of steepest 
descent (MSD) 0] (saddle-point approximation) for per- 
forming an occurring contour integral, and (ii) the as- 
sumption of a small kinetic momentum of the escaping 
electron. The Keldysh approximation and variants of 
it are also very popular, because they are the basis for 
methods that predict strong-field ionization rates also for 
heavier atoms and molecules (see, e.g., [3, Q). 

Recently, there has been proposals to obtain an exact 
Keldysh theory by removing the two additional simplifi- 
cations (i) and (ii). It was especially suggested to avoid 
approximation (i) by solving the occurring contour in- 
tegrals with the aid of the exact residue theorem (RT) 
instead of the MSD @, 0, S]- It was shown that the RT 
method yields a two times larger transition amplitude 
for the Is state of a hydrogen-like atom and, as a conse- 
quence, a four times larger ionization rate. Furthermore, 
the ionization rate of the first excited (2s) state obtained 
in [f| when applying the RT differs significantly from the 
MSD result. In view of the popularity of the MSD ap- 
proximation for treating strong-field problems like ioniza- 
tion Q or high-harmonic generation [l(| this is of course 



a very important result. This has motivated the present 
study in which a careful reinvestigation of the RT and 
the MSD is performed (Sec.|TI]). It is shown that the ap- 
plication of the RT as proposed in 0, 0, Q contains an 
unjustified neglect of the contribution of one integral and 
that it is this omission which is the main source for the 
previously reported deviation between the RT and the 
MSD results. Therefore, the MSD provides in fact more 
reliable results than the (incomplete) RT approach. It 
is furthermore discussed that the MSD fails for Rydberg 
states and a correction to it is proposed. The conclusions 
of Sec. |IT]are supported with the aid of a numerical study 
in Sec. EH 



II. THEORY 

A. Transition amplitude 

In order to provide the basis for the subsequent discus- 
sion and to introduce the notation a brief summary of the 
Keldysh theory is given that follows closely the one de- 
scribed in Appendix A of The total ionization rate 
of a one-electron atomic system with the electron bind- 
ing energy in the harmonic laser field F(i) = F coswi, 
with the period T = 2n/u!, can be expressed as the sum 
over Af-photon processes (atomic units are used through- 
out this work) 



W = 2ir 



d 3 p 
(2^p 



|A(p)| 2 £ 5{E h + P -+U p -Nuj) (1) 



N=N„ 



where U p = F 2 /(4a; 2 ) is the electron quiver (ponderomo- 
tive) energy due to the field. The transition amplitude 
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A(p) can be calculated using 



i 

A(p) = ^J dt J d 3 r%(r,t)V F (t)H> (T,t) (2) 



where ^o(r,t) = e iBbt $o( r ) is the wave function describ- 
ing the initial electronic state in the atomic potential 
U(r). Therefore, $o( r ) fulfills the stationary Schrodinger 
equation 



iv 2 + C/(r) 



*o(r) = 



(3) 



The interaction with the laser field is given in length 
gauge by Vp(t) = r ■ F(i). Finally, the Volkov wave func- 
tion [12j \& p (r, i) satisfies 



i 



v 2 + V F {t) 



and can be explicitly written as 
* p (r,i) = exp 



ivrz{t) - / 7V 2 {t')dt' 



(4) 



(5) 



where 7r(t) = p+ (F/w) sinwi is the mechanical momen- 
tum of an electron with the canonical momentum p in 
the field F(t). Introducing the auxiliary functions 



^o(q) 



d 3 re~^ r (F ■ r)4> (r) = *F • V q <l (q) (6) 



(where <I>o(q) is the Fourier transform of &o{r)) and 
t 



S(t) = J dt' 



(7) 



the transition amplitude A(p) can be rewritten as 

T 

A{p) = y I " dt C0S M v ° (*(*)) etS(t) ■ 



(8) 



The equivalence of A(p) in and L(p) in Eq. (15) 
of the original Keldysh work [2| can be shown in the 
following way. The path of the integration over i in (J5J) 
can be shifted into the complex plane by means of the 
transformation t = t + it where e is an infinitesimally 
small positive number. Introduction of the new complex 
variable u = sin Lut = sin tot + it cos tot transforms the 

T 

integration J dt to one on the closed contour C; n which 
o 

encloses the interval (—1,1) (see Fig.[T]a). Applying the 
same procedure to the integral contained in the function 
S(t) yields 



A{p) = j duV (u)e lS{u) 
a, 



(9) 
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FIG. 1: (a) The path Cm of the contour integration for A(p). 
(b) The path C u of the contour integration for S{u). (c) 
The asymptotic behavior of exp[iS(u)], contours C± around 
saddle points w±, contour C ou t, and deformed Ci n used within 
MSD. 
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V (u) 
S(u) 



7R 



#b 1 / F 
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to Zlu \ to 



(10) 
(11) 



The transition amplitude A(p) in ^ is for f(v) — 
Vl — t> 2 identical to £(p) in 2] which is also the start- 
ing point of the analysis in [6]. Since the square root is 
usually defined to possess a non-negative real part which 
can in the present context be misleading and cause a sign 
error, we introduce the function 



f(v) = Sign[Im(v)] yl — v 2 , 



(12) 



which is analytical in the whole complex plane except its 
branch cut [—1, 1]. The path of integration C u specifies 
the clockwise path around the branch cut (see Fig.fTJb) 
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starting at v = ie and terminating at v = u. Note, S(u) 
is a multivalued function, so we have selected also the 
branch cut along positive imaginary axis. 

Due to the delta function in Eq. (Q} one needs to calcu- 
late ^(p) only for |p| = p N = y/2(Nuj - E h - U p ). Both 
© and 1(5)) can equivalently be used for numerical inte- 
gration to yield an exact result. The use of (JSJ) provides 
more flexibility, since the contour C- m can be deformed in 
a convenient way. 

There exist two special points u± (|A4p in the complex 
plane u. They are simultaneously the saddle points of 
S{u) and poles of V (u). In @ Keldysh has used MSD 
to approximate A(p). In Q the authors have proposed 
to use RT Q for an exact calculation of ^.(p). For the Is 
state of hydrogen-like atoms the expression for A(p) ob- 
tained in [6( is larger than that of Keldysh using MSD by 
exactly a factor of two, provided the small p approxima- 
tion is consistently used or omitted in both the RT and 
the MSD approach. As is shown below, the disagreement 
is a consequence of a wrong assumption made in [6J. In 
fact, for excited states of hydrogen-like atoms the therein 
proposed approach may lead to drastically wrong results. 

For spherically symmetric bound states of hydrogen- 
like atoms with principal quantum number n the function 
Vq(u) can be presented (see Appendix B) as 



c' s 



V (u) = 



9+{u) = g-{u) 
(u — u+) v (u — u-) 1 



(13) 



where g±(u) = Vq[u)[u — u±) u is an analytical (and, 
possibly, slowly varying) function in the vicinity of the 
special points u±. Both procedures (MSD and RT) are 
considered in this work for general Vq (u) having poles of 
order v at u = u±. 



ct 



C 



FIG. 2: The integration paths Cf 1 passing around the special 
point u s from both sides and given parametrically in Eq. (|20|l . 
The contours CI and C\ are used to connect the contours 
at infinity. 



In [6[ the value of I out is implicitly assumed to be zero. 
However, a simple analysis shows that there are no rea- 
sons for such an assumption. Indeed, for u — Re with 
R — ► oo, one finds [see Eq. (|A2|) ] 

e *S(u) _> e -(U p /uj)R 2 COS {26) e -(U p /u)R 2 S \n{2e)i^ /jgN 

Since Vq{u) ~ i? -5 for R — > oo, the integrand in (jTBJ) has 
the following limits: 



\V {u)e iS ^\ 
\%{u)e iS ^\ 



7T 7T 3-7T 57T 

°' -A <6< r- <d< -' 

3lT 7T 7T „ 37T 

°°! — r < 9 < t < 6 < ~r- 

4 4 4 4 



Therefore, it is impossible to select C ou t in such a way 
that the integrand on the whole contour approaches zero. 
Moreover, in section IIIII it is numerically demonstrated 
that I ou t is of the same order of magnitude as A(p) or 
even larger. 



B. The residue theorem 

Since for |p| = the function Vo(u) exp[iS(u)} is ana- 
lytical in the whole complex plane except the branch cut 
[—1,1] and the poles u±, © can be modified using 



(14) 



C in C + C- Cout 



where C± are contours around u± and C ou t is a contour 
enclosing (—1, 1) and u± (see Fig.[T]c). The integrals A± 
along C± can be calculated using RT which yields 



-4. 



2™ 



lim 



v-l 



5±(«)e lS(u) 



(y — 1)! u^u± du u 1 
With the knowledge of the integral along C ou t , 



duV Q {u) e lS(u \ 



c out 



the value of A(p) can be calculated using (fT4| as 
A{p) = (A + + A_) - 7 out . 



(15) 



(16) 



(17) 



C. Contours through steepest descent 

In order to understand the appearance of the factor 
two between the RT and the MSD results it is important 
to connect the two approaches. For this purpose, the four 
auxiliary integrals 



9s(u) 



(«-«»)" 



-du 



(19) 



are introduced where the index s = ± specifies one of the 
two special points. The contours Cf- are given paramet- 
rically by 



u s + (x ± ie) 



-oo < x < oo, e — >• 0+ (20) 



starting at x — > — oo. Here, plus (minus) corresponds 
to the contour passing slightly above (below) the spe- 
cial point u s (see Fig. ^ . The value of Q s is chosen in 
such a way that the contours are passing through the 
steepest descent of iS(u), i.e. as 



2i 



S"{u s ) 



(21) 
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where the argument Q s satisfies — it /A < argQ s < 7r/4 
(see Eas.(|X8jl). According to |TSJ the integrand in (JTOJ) 
must then exponentially decay to for x — > ±00. This 
specific choice of Q s allows to directly apply MSD in the 
following subsection. Here, Cauchy integration rules Q 
are employed to deduce three useful relations: 

(i) Deforming the contour in A(p) to pass along Ct 
in positive direction, along C^T in negative direc- 
tion and connecting the ends of these contours at 
infinity one obtains 



A(p) =17 -It 



(22) 



(ii) Connecting contours Cf at infinity with the con- 
tours C l s and CJ (see Fig. [2]) and applying the 
residue theorem one obtains 



A s = IT - It 



(23) 



(iii) Substituting and @ into (T7]) one obtains 



lout — I- — I-L 



(24) 



Note, the use of a sufficiently small but finite positive e 
yields the same value of if. This is used to compute if 
numerically 

Keeping in mind that equations (|22p, (f2"3"|) , and (j2~4"|) 

are exact and no approximations have been done so far, 
we apply now MSD to approximate if. 



D. The method of steepest descent in the presence 
of a singularity 

Since no difference is made between two different con- 
tour integrations around the same saddle point (as, e. g., 
for /+ and I~) in [11] (Appendix B), we shortly repeat 
the main steps. From (|20p the relations 



du = Q s dx, 



1 



(±iy 1 



{u x -u s y (iQ s y(tTix) v 



(25) 



follow. We expect the vicinity of u s to give the main 
contribution to the integral and assume that g s {u) is a 
slowly varying function in the vicinity of u s . Then, using 
the approximation 



g s (u)e lSM « g s (u s )e iS( - u >h- x2 * 2iex+e 
and the identity 

00 

1 1 



(26) 



(a Tib)" T(u) 



drir) v - 1 e- r > a e ±ir > b , a > (27) 



one obtains 



(±1) 



g s (u s )e iS ^ 



v-X -r)e+e 



drj rf e 



00 

x J dx e 

—00 



dxe- x2±l{v - 2t)x . 



(28) 



The integration over x and 77 yields 

irg s {u s y s( - u ° 



(±1)"J„ 



i"Qr x r(^i)' 



(29) 



Therefore, MSD predicts if to be equal for even v and 
to differ only by the sign for odd v. Using equations (|22|) 
and (12411 this result can be rewritten as 



/out = {-\) v ~ x A{p) [within MSD]. (30) 

Using ([221) the prediction of MSD for A(p) is 

A MS D = (-iyi+-I- ■ (31) 

Substitution of ([30]) into (T7]) shows that for odd v MSD 
yields 

A(p) w (A+ + A_)/2 [MSD, odd 1/]. (32) 

Its value is thus two times smaller than the one obtained 
with the assumption / out = 0. For even v MSD predicts 



|A(p)| > \(A+ + A-)\ [MSD, even v\ 



(33) 



Note, ([3^|) and (|3"3"]) are valid for every Vq(u) satisfying 
(|13p. if g±(u) is a slowly varying function in the vicinity of 
u±. The fact that for the Is state of hydrogen-like atoms 
[y = 3) one finds exactly a factor 2 difference between 
MSD and RT and thus an equality sign in (|3"2")l should 
be seen as an accidental case that is due to the relative 
simplicity of Vq(u) for the Is state. 

In section lnil MSD is tested numerically and it is shown 
that the assumption of a slowly varying function g±(u) 
is valid only for small n (or v). 



III. NUMERICAL TEST 

To support our conclusions of the previous section the 
example results of a numerical study are reported in 
Table HI For a more transparent analysis a number of 
parameters were fixed. This includes the amplitude of 
the electric field F = 0.02 a.u., the frequency u> = 0.01 
a.u., the binding energy E Q — Z 2 /(2n 2 ) = 0.5 a.u. (thus 
charge Z — n), the angle F • p = 0.9, and the number 
of photons N — 161. With such a choice of fixed pa- 
rameters a variation of the principal quantum number n 
leaves the function S(u) unchanged (see Appendix A for 
details). The same is true for the positions of the special 
points u± and the values of Q±. Therefore, only function 
Vo(u) varies with n. 

We use (j8]) to calculate the exact value A cx of A(p) 
and (fT5|) to evaluate A± . The integrals if are calculated 
numerically, Eq. (|22|) and (|23[) are used for a check of the 
numerics and 7 ou t is obtained from (|24D . 

As follows from the discussion above, the condition 
I~ w (— l) u If obtained in (|2"9")) can be used as a cri- 
terion for the validity of the simple MSD formula (|3ip . 
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TABLE I: Contour integrals [Eq.JT9j], J ou t [Ea.[f2l|)]. quantities I ± [Eq.(T29j)], and A+ + A- [Eq.JT5j] for different 

principal quantum numbers n and fixed parameters (F — 0.02 a. u., ui — 0.01 a. u., Eh — 0.5 a. u., F ■ p = 0.9, TV = 161). The 
exact value A ex for the amplitude A(p) [Eq,(|8jl] is compared with the prediction of the simple MSD formula Amsd [Eq. (l31|) ]. 
the corrected MSD formula A c msd [Eq. (IC5|l ]. and the amplitude Akm (with appropriate phase normalization) given by the 
"two-term saddle-point approximation" |13j ]. 





n — 1, x 10 8 


n = 2, xlO 7 


n = 3, xlO 6 


n = 4, xlO 5 


n = 5, x 10 


n 


-0.212 - 1.560 i 
0.226 + 1.883 i 
-0.219 - 1.722 i 


-0.320 - 1.760 i 
-0.302 - 1.884 i 
-0.312 - 1.812 i 


-0.374 - 1.652 i 
0.274 + 1.347 i 
-0.325 - 1.495 z 


-0.371 - 1.370 i 
-0.189 - 0.758? 
-0.276 - 1.045 i 


-0.326 - 1.034 i 
0.105 + 0.352? 
-0.201 -0.644 i 


ii 

IZ 

I- 


-1.180 - 1.485 1 
0.960 + 1.248 i 
-1.070 - 1.366 i 


1.127 + 1.539 i 
1.027 + 1.465 i 
1.069+ 1.496 i 


-0.765 - 1.142 1 
0.912 + 1.428 i 
-0.835 - 1.282 i 


0.406 + 0.667 i 
0.713 + 1.227 i 
0.549 + 0.931 i 


-0.176 - 0.322 i 
0.506 + 0.959 i 
-0.317 - 0.595 i 


A-+A+ 

lout 


2.578 + 6.175 i 
1.172 + 2.808 i 


-0.083 - 0.199 i 
1.346 + 3.224 i 


2.325 + 5.569 i 
1.286 + 3.080 i 


0.490 + 1. 172 i 
1.084 + 2.597 i 


1.113 + 2.666 i 
0.832 + 1.993 i 


Amsd 

A C MSD 

A ex 
Akm 


1.289 + 3.088 i 
1.406+ 3.367 i 
1.406+ 3.368 i 
1.378 + 3.300 i 


-1.381 - 3.308 i 
-1.432 - 3.431 i 
-1.429 - 3.423 i 
-1.680 - 4.025 i 


1.160 + 2.778 z 
1.047 + 2.506 i 
1.039 + 2.489 i 
1.641 + 3.932 % 


-0.825 - 1.976 i 
-0.603 - 1.444 i 
-0.595 - 1.425 i 
-1.365 - 3.270 i 


0.517+ 1.239 i 
0.287 + 0.687? 
0.281 +0.673 i 
0.999 + 2.393i 
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FIG. 3: (Color online) Demonstration of the failure of the 
MSD for large principal quantum numbers n. The ratio 
|/out/A cx | (boxes) grows exponentially with n which causes 
an increase of the ratio IAmsd/AexI (circles) between the ap- 
proximate and exact amplitudes Amsd and A ex , respectively. 
Therefore, the simple (standard) MSD formula [Eq. (|31|l ] fails 
and must be improved to be applicable for large n. The in this 
work proposed corrected MSD formula [Eq. (|C5|l ] shows very 
good accuracy in a large region of n, as can be seen from the 
corresponding ratio |A c msd/A gx | (triangles). The parameters 
used in the computation are the same as in Table U 



Table U shows that this condition is fulfilled for n = 1, 2 
and the relations (|30I32I33[) are valid. With increasing n 
the condition is, however, not well fulfilled and the accu- 
racy of the MSD prediction decreases. Numerical tests 
show that for large n the values |/+|, \IZ\ are by orders 
of magnitude larger than \I+ |, \I^\. This leads to the 



following relations (see Fig. [3]) 

|A CX | < |A MS d|, Kx| < |4ut|, lout ~ A+ + A— 

Therefore, function Vq(u) cannot be given as simple as 
in (fl3|) . Instead, Vq(u) can be represented in the vicin- 
ity of u s by a sum over terms having different orders of 
poles (see Eq. (|B7|0 . Moreover, it is possible to consider 
also higher derivatives of S(u), as is done in Eq. (|C2p . 
The resulting representation of Vq(u) exp[iS(u)] given in 
Eq. (|C3|) and the subsequent use of MSD to it leads to 
a much higher accuracy. As can be seen from Table U 
and Fig. [3] the corrected MSD formula (|C5[) yields a sig- 
nificant improvement and can be used for the numerical 
computation of A(p) in a large range of n. This paves the 
way for a detailed study of the validity of the Keldysh ap- 
proximation for, e. g., Rydberg atoms that is not blurred 
by a failure of the usually adopted MSD approximation. 

It is instructive to compare A c msd with the rece ntly 
published "two- term saddle-point approximation" [13j. 
This approximation is a modification of the "one-term 
saddle-point approximation" [TlT |. Both approximations 
are based on the principle that only the leading term of 
the Laurent expansion of the Fourier transform $0(0.) 
at the saddle points is considered. Since in the "one- 
term" approximation the contributions from higher-order 
derivatives of S(u) are also ignored, it is equivalent to 
the simple MSD formula Amsd and yields thus the same 
numerical results. (Note, the resulting expressions are, 
however, different, because different complex variables 
(f> = u>t and u are used in the derivations.) As can be 
seen from Table HI the amplitude Akm yielded by the 
"two-term" approximation (using an appropriate phase 
normalization) gives relatively good agreement for n = 1, 
although it is clearly less accurate than A c msd- For the 
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excited states the "two-term" formula fails rapidly and 
yields results which are even less accurate than the ones 
obtained with the "one-term" formula. To understand 
this fact, we remind that the second term in [l3[ comes 
from the next-to-leading pole, and only the contribution 
arising from the third derivative of S(u) is considered. 
A detailed analysis shows that a singularity of the same 
order can also arise from the next-to-leading term of the 
Laurent expansion of $o(*l)- I n addition, for larger values 
of n it can also be important to consider further terms of 
the Laurent expansion. Since the corrected MSD formula 
(|C5[) takes all of this into account, the resulting A c msd 
is significantly more accurate than Akm- 



APPENDIX A: CALCULATION OF u±, 

S(u),S'(u),S"(u±),S"'(u±), AND Q±. 

The integration of (fTl"]) for |p| = pn yields 



exp[iS(u)] = exp 



pF r 



.[l-/(u)]-i-Eu/(u) 



x [fiu) + iu 



N 



(Al) 



For u = Re with R — > oo, one has f(u) — > —iu + iu 1 /2 
and 



exp[iS(u)] 



(-2iu) 



N 



(A2) 



IV. CONCLUSION 

In this work it has been demonstrated that the residue 
theorem was not correctly employed in [y, LZ| , since the 
derivation was based on an unjustified assumption that 
one integral vanishes. This neglected term is, however, 
of the same order of magnitude as the remaining ones or 
even much larger. For the Is state of hydrogen-like atoms 
it is almost identical, and thus its omission results in an 
overestimation of the transition amplitude by a factor 
two for this case. This deviation was in t (| incorrectly 
assumed to be a failure of the widely used saddle-point 
approximation. Considering a 2s state, it is furthermore 
concluded that an application of the method proposed 
in 0] to a 2s state would yield an even larger (erroneous) 
deviation. 

Such a large deviation for the 2s state was in fact re- 
ported in Q where also the residue method had been 
applied. Analogously to Q the deviation was attributed 
to a failure of the saddle-point method, but is in fact 
due to the same unjustified omission of a non- vanishing 
integral. The direct applicability of the present findings 
to [H can be verified, since the derivation in Q differs 
from Q essentially only by the choice of 4> = ut as com- 
plex variable, while in [6[ and the present work u — sin uit 
was used. 

The applicability of the method of steepest descent 
(saddle-point approximation) for arbitrary ns states has 
also been investigated in the present work. It is found 
that the simple standard formula fails for large n. To 
overcome this problem a corrected formula is proposed. 



The saddle points u± of S(u) can be determined by 
the following condition: 



E h If F , 
— + —[p+-u =0. 

uj 2u> \ uj 



(A3) 



Introducing the Keldysh parameter 7 = klo/F with k = 
V2i/b, the scaled momentum x = Pn/k, and £ = F • p 

u± =-a±pi, a = 7xC, P = 7V 7 ! + X 2 (l - C 2 ) (A4) 

is obtained. Using (|A4p the first derivative S'(u) can be 
expressed as 



S'(u) = 



2U P (u — u + )(u — «_) 
uj f(u) 



(A5) 



and the values of the second S"(u) and third S'"(u) 
derivatives at u — u± are given by 

b [u±)=±— — r, S (u±) = — ipr- — . (A6) 



"w/(w d 



ujf 3 {u±) 



The absolute value Q and the argument of Q s de- 
fined by (f2Tj) can be written as 



Q = 



2cj7 



[(l + a 2 +p 2 ) 2 -4a 2ll/ 



(A7) 



tan4?9 ± = ± 



2ap 



l + p 2 -cr 2 ' 4 " " ± ~ 4 



7T 7T 

< ^+ < -. (A8) 
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APPENDIX B: FUNCTION V (u) FOR THE nS 
STATES OF A HYDROGEN-LIKE ATOM. 



Then for g± (u) denned in ([13 



Consider the spherically symmetric state (with prin- 
cipal quantum number n) of a hydrogen-like atom with 
potential Z7(r) = Z/r, where Z is the charge of the nu- 
cleus. Its Fourier transform is given by 

*o(q) - -A £(-l) fc 2 2 *C«+* (-Ji-z) (Bl) 



fc=0 



where are binomial coefficients. Using the identity 
^q/(<7 2 ) = % c ldf(q 2 )/(dq 2 ) one can rewrite Eq. © as 



3±H= J2 M - 



(n,n+2— m) 



(u-u±) w , (BIO) 



m=Q 



one has 



9± (u ± )=Mt n ^ = (±2r-^ in+lh2n+3F 



(Bll) 



APPENDIX C: CORRECTED MSD FORMULA 
FOR THE nS STATES OF A HYDROGEN-LIKE 
ATOM. 



where 



^ = (-l) fc (fc-l)2 2fc - 4 C"+_ fe 5 - 3 



Introducing 



and using 



P±(«) 



7 



u — u± 



(B3) 



(B4) 



F q = ^[P-^^ + PI^u)}, = P+(u)P_(u), 

2 q z + k z 



Representing exp[iS'(u)] as 

exp[iS(u)] = exp[iS'(u) - (i/2)S"'K)(u-M s ) 2 ] 

x exp[-g 2 (u- Ms ) 2 ] (CI) 

and performing a Taylor expansion of the first term on 
the right hand side of (|C1|) at u — u s yields 



e iS(u) _ e iS{u 3 ) e -Q 2 3 {u-u 3 ) 2 



x 1 



iS"'{u a ) 



{U-U s f + 



(C2) 



Keeping the first two terms of the Taylor expansion the 
integrand in (|9I16[) can be rewritten as 



one can rewrite Eq. (|10[) as 

n+2 



V (u) = J2 D^{Pl~ 1 P , l + P^- 1 }. (B5) 



fc=3 



Introducing R± = P±(u T ) = ^fj/(2pi) and using the 
Taylor expansion of P± (u) at u ~ u T , 

oo 

P|(u) = £(-l) m C^ +m - 1 i?^ +m p-">), (B6) 
one can rewrite Vq (u) as a Laurent series at u ~ u± 



n+2 M^ n,z ^ 



U— — OC 



{u-u±) u 



(B7) 



where 



n+2 



r— max(f ,3) 



, V — — oo 



(it - U s ) 1 



6 ^ (u - u s ) v 



(C3) 



Omitting terms with negative ^ in (|C3[) and applying the 
procedure described in Sec. Ill Dl 



a?'"(u.) ^ (±l)»M<""' +3 > 1 

is obtained as approximation for 1^ . Note, that neglect- 
ing all terms with v < n + 2 in (|C4p and using (|B11|) 
one obtains the simple (standard) MSD formula (|2U)) for 
Jg . A corrected approximation for A(p) is then obtained 
using (|2"2"|) as 



and 



Q« = (-l) m {C* +m - 2 - C£t?- 2 }, Q { k) = 1. (B9) 



— I+-CMSD ~ I-;cMSD- 



(C5) 
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